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Abstract 

Nonlinear plasma oscillations in an arbitrary mass ratio cold plasma have 
been studied using 1-D particle-in-cell simulation. In contrast to earlier work 
for infinitely massive ion plasmas it has been found that the oscillations phase 
mix away at any amplitude and that the rate at which phase mixing occurs, 
depends on the mass ratio ( A = m_/m + ) and the amplitude. A perturba- 
tion theoretic calculation carried upto third order predicts that the normalized 
phase mixing time uj p -t m i x depends on the amplitude A and the mass ratio A 
as ~ [(A 2 /24)(A/\/l + A)] -1 / 3 . We have confirmed this scaling in our simu- 
lations and conclude that stable non-linear oscillations which never phase mix, 
exist only for the ideal case with A = 0.0 and A < 0.5. These cold plasma 
results may have direct relevance to recent experiments on superintense laser 
beam plasma interactions with applications to particle acceleration, fast ign- 
itor concept etc. 

PACS number (s): 52.35.Mw, 52.65.Rr 
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The physics of the damping of nonlinear cold plasma oscillations is a topic of considerable 
fundamental interest since it is the simplest nonlinear collective irreversible phenomenon 
characterising the plasma state. It also has wide applications to a number of problems 
of current interest such as particle acceleration by wakefields and beat waves created by 
intense lasers or particle beams, the fast ignitor concept in inertial fusion where relativisti- 
cally intense coupled electromagnetic - plasma waves modes propagate deep into overdense 
plasmas to create a 'hot spark' and a number of other astrophysical / laboratory / device 
based plasma experiments where intense plasma oscillations are generated. The conventional 
thinking about the physics of this interaction is well illustrated by the exact solution for 
nonlinear one - dimensional cold plasma fluid equations with infinitely massive ions. These 
exact solutions may be obtained by transforming to Lagrangian coordinates as shown in 
]l|-[3|] or using stream functions [|J. The exact solution shows that coherent oscillations at 
the plasma frequency u p are maintained indefinitely over the region of initial excitation, 
provided the normalized amplitude of the initial density perturbation A(= — ) is kept be- 
low 0.5. For A > 0.5, one expects and observes wave breaking and fine scale mixing of 
various parts of the oscillation ]I] . Mathematically, the electron number density blows up 
at A = 0.5; this is because the Jacobian of transformation from Eulerian to Lagrangian 
coordinates goes to zero as A — > 0.5 and the transformation is no longer unique. Physically, 
this is equivalent to crossing of electron trajectories which leads to multistream motions and 
wave breaking as discussed in . Studies of wave breaking and phase mixing damping are 
based on numerical simulations. 

The above description is adequate when the background positive species are infinitely 
massive = A — ► 0) and are uniformly distributed. If the background is inhomogeneous, 
then as was shown by Dawson |L|, cold plasma oscillations phase mix away in a time scale 
t ~ 2(Aj 7 }dx)x ' a ^ arbitrarily low amplitudes. For a sinusoidal distribution of background 
species, such a phenomenon in the form of mode coupling of a long wavelength mode to 
short wavelength modes was observed by Kaw et. al. ||. They found that, the time scale in 
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which energy goes from long wavelength mode to short wavelength mode is t ~ — — , where 
V is the amplitude of the background inhomogeneity. The exact solution for the cold plasma 
oscillations in a fixed sinusoidal background was given by Infeld et. al. who described 
phase mixing in terms of electron density burst. 

In this paper we show that the phenomenon of phase mixing will also occur in a homo- 
geneous plasma at arbitrarily low amplitudes, provided the background positive species are 
allowed to move (A / 0). This is because the background species respond to ponderomo- 
tive forces either directly or through low frequency self-consistent fields and thereby acquire 
inhomogenities in space. Such an effect has been observed in electron positron plasmas 
(A = 1) by Stewart [0]. In plasmas with finite temperature, it is well known that plasma 
waves dig cavities by ponderomotive forces and get trapped in them; this is the physics of 
strong turbulence of Langmuir waves as elucidated by Zakharov and leads to envelope 
soliton formation in one - dimension and collapse phenomenon in 2-D and 3-D. In a cold 
plasma, stationary states cannot form even in 1-D because there is no thermal pressure effect 
to counterbalance the ponderomotive forces. The result is that the density cavities being 
dug by plasma oscillations have an amplitude which increases secularly in time . Similarly, 
the response of plasma oscillations to the presence of density cavities is also different from 
that of Zakharov problem. In the Zakharov problem, the thermally dispersive plasma waves 
get trapped in density cavities forming localized wavepackets. Here, the inhomogeneity of 
the cold plasma ( due to the self - consistently generated perturbation ) causes different 
parts of the plasma oscillation to oscillate at different frequencies || ^ @ resulting in in- 
tense phase mixing of plasma oscillations. Thus we physically expect that if the background 
species is allowed to move and get redistributed into inhomogeneous clumps of density, the 
phase mixing damping of cold plasma oscillations should come in at any amplitude and is 
not restricted to waves with A > 0.5. It may be emphasised here that for many applications 
involving interaction of superintense laser beams with plasmas ( such as particle accleration 
by wake fields, penetration into overdense plasmas etc. ) the cold plasma limit considered 
by us is more relevant than the Zakharov description, because typically the plasma wave 
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intensities are such that | E \ 2 /AnnT 3> 1. 

In this paper we carry out particle simulations for elucidating the physics of phase mixing 
damping of nonlinear cold plasma oscillations in an arbitrary mass ratio plasma (A arbi- 
trary). We also present a perturbation - theoretic analysis to give a quantitative estimate of 
the phase mixing time for moderate amplitude oscillations and compare it with simulation. 

We start with the cold plasma equations viz. the continuity equations and the equations 
of motion for the two species and the Poisson equation. We introduce new variables V, v, 5rid 
and Sn s defined as V = v + + 1>_, v — v + — t>_, 5rid = Sn + — 5n^ = n + — n_ and 5n s = 
5n + + 5n_ = n + + n_ — 2 to write the cold plasma equations in the form 

9 lSnd + dllv+ Y^i±l^ ] = (1) 

aj n . + a,\y+ VSn ' + vSn ' ] = o (2) 

V 2 + V 2 

d t V + d x {^—) = -(1 - A)E (3) 

d t v + d x (^) = (l + A)E (4) 

d x E = 5n d (5) 

Note that we have used the normalizations: n± — > n±/no, x — > kx, t — > u p ^.t, v± — > 
v±/uj p -k~ 1 , E — > E / (AnriQek" 1 ) , with uj 2 _ = 4nn e 2 / m_ and A = m_/m + . 

Using n-(x, 0) = 1 + 5 cos /ex, 0) = 1 and i>±(x, 0) = 0, as initial conditions, the 
solutions of the linearised equations are 

5n d ^ = Acoskxcosujpt (6) 

E^ = — sin kx cos u p t (7) 

drSp = - — — A cos — cos cjpt) — A cos fcx (8) 

y(i) = - 1 A A sin £;:r sin w p t (9) 
kujp 

v (i) _ ^ ^4 sin fex sin a> p t (10) 

where A = —5 and cj 2 = 1 + A. At this level of approximation, the solutions show 
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coherent oscillations at the plasma frequency u p . Both the species oscillate with the same 
frequency which is independent of position. 

In the 2nd order, the solutions are expressed as: 

5n^P = — A 2 cos 2kx[- — (— + -u p t sin u p t + - cos2uj p t — cosu p t) — -u p tsmuj p t} (11) 

1 | / * ^4 — 4 

A 2 1 — A 1 1 1 1 

SE^ 2 ' = -sm2kx[ — ( — h -u p tsmu p t H — cos2u;pt — cosuj p t) — -uj p tsm.u p t\ (12) 

2/c 1 ~\~ 2 4 2 4 

. (2) A 2 , . A 2 A(l - A) 3, 

^ = -y cos2fcx[ 1 - - -j—-^^u p t sm Upt- -(1- cos 2u p t) - (13) 

.l-A.o, 5 11.. 

(Y+~K' (2 cos u p t - - cos 2u p t - —)] 

,r(o\ A 2 , r A w-.l-A.,, 1-3A 

^ = __ sm2fa [_ t + ^(—^(jjj-^oo.^ + (14) 

7-5A . 5 . n , 1 . „ , 
■ sm Upt — - sm 2uo p t) u p sin 2u p t\ 
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2(1 - A) 4 " 8 

^ ~ rr^ sin2fcx[sinu;pt — WptcostUpt — —(2sm2uj p t — 3smuj p t — uj p tcosuj p t)] (15) 

o/c 1 Za 

The 2nd order solutions clearly exhibit generation of 2nd harmonic in space and time 
as well as bunching of plasma particles in space. Both of these features are also evident in 
the solution of Kaw et. al. || and Infeld et. al. ||; but in contrast to their work, where 
the background ion density was kept fixed in time, here the density of the plasma particles 
self - consistently changes with time as ~ t 2 , as seen in the expression for 5nf\ Because 
of variation of plasma density with time, the phase mixing of an initial coherent oscillation 
happens much faster in this case. To make an estimate of the phase mixing time, consider 
the charge density equation ( 5n d in this case ). The equation for 5n d correct upto third 
order stands as 

d tt 5n d + u 2 p [l + ^(5n« + 5nW)]5n d « (16) 

In the above equation, if we neglect the 2nd order term, then we essentially get the same 
phase mixing time as in Ref. || modified by a factor which depends on A. Now taking only 



the leading order secular terms from the expressions of dny' and 5nf^ ( there are no secular 
terms in 5n^ ) we get 



d tt 8n d + ujl[l + 



A 2 t 2 A 

Aoj 2 
p 



cos2kx]5n d ps 



(17) 



Using the initial conditions 8n d = A cos kx and d t 5n d = the WKB solution of the above 



The above expression clearly shows that the energy which was initially in the primary 
wave at mode k goes into higher and higher harmonics as time progresses. This can be 
interpreted as damping of the primary wave due to mode coupling to higher and higher 
modes. Microscopically, as the plasma particles oscillate at the local plasma frequency, they 
gradually go out of phase and eventually the initial coherence is lost. Because of generation 
of higher and higher harmonics with time, the charge density becomes more and more spiky 
and as a result the electric field gradients become more and more steep. This does not go 
on indefinitely. In reality, the density peaks get limited by thermal effects with the Landau 
damping of high k modes by resonant particles coming into picture. This process takes 
energy from the high k modes and puts it on the particles, thereby raising their temperature, 
which in turn limits the density peaks by exerting a pressure gradient. The time scale in 
which the initial coherence is lost ( or the phase mixing time ) can be seen from equation 
(18) as uj p -t m i X scale as ~ L4 2 A/(24a/1 + A)] -1 / 3 . It shows that only for the ideal case 
A = 0.0 (infinitely massive ions), phase mixing time is infinity, i.e. the initial coherence is 
maintained indefinitely 0|| . For an actual electron - ion plasma, A although small, is finite 
and hence plasma oscillations in it phase mix away at arbitrarily small amplitudes and in a 
time scale dictated by the amplitude of the initial perturbation. 

Now we present results from a 1-D particle-in-cell simulation which confirms our scal- 
ing of phase mixing time. For numerical simulation, we have used a 1-dimensional model 
with periodic boundary conditions and have followed 5120 electrons and as many positively 



equation is 




(18) 
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charged particles ( the plasma taken as a whole is neutral ) in their own self - consistent 
fields. The particles are initially at rest and the system is set into motion by giving a density 
perturbation of the form n_ = 1 + 5 cos kx to the electrons. In the simulation, we follow the 
time development of various modes of charge density (Snd) . To compare with our theoretical 
model, we rewrite equation (18) as 

A n=oc nir 
5rid = — J n [a{t)][cos(ujpt + — ){cos(2n + l)kx + cos(2n — l)kx} + (19) 

n=— oo 

T17T 

sm(uj p t + — ){sin(2n + l)kx + sin(2n — l)fcx}] 

where a(t) = 4^--^==. The amplitude of the first fourier mode can be seen from the 
above equation as 

| 5n d | n=1 = |[J 2 (a(t)) + J?(a(t))]3 (20) 

It is clear from equation (19), that upto the order of approximation considered, there 
are no even number modes in the system. Fig. (1) show temporal variation of | 5nd | n =i for 
A = 1.0 and A = 0.05. The dotted curve is the simulation result and the solid line shows 
our expression (20) for the envelope of the oscillations. It is clear from the figure that our 
approximate expression (20) captures the early evolution of the plasma quite well. Fig. (2) 
shows the variation of T m i x = uj p t m i X with A for a fixed A = 0.1 (curve(l)) and with A for a 
fixed A = 0.01 (curve(2)). These curves clearly confirm our formula for phase mixing time. 

In conclusion, we have demonstrated that nonlinear plasma oscillations in a cold ho- 
mogeneous plasma, phase mix away at arbitrarily low amplitudes. This is because during 
the course of motion the plasma particles respond to ponderomotive forces, acquiring in- 
homogeneity and thereby making the plasma frequency a function of space. As a result, 
electrons at different locations oscillate with different ( local ) plasma frequencies and the 
imposed plasma wave losses coherence. The formation of density clumps can also be seen 
from the Zakharov equations || for a warm electron - ion plasma. According to Zakharov, 
the slow variation ( in the ion time scale ) of the background density in the presence of a 
high frequency oscillation is governed by duSn s — Td xx 5n s = d xx \ E \ 2 . In the limit when 
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the thermal term balances the ponderomotive force term ( i.e. ^ ~ ), we get caviton 

solutions in 1-D which are unstable to transverse perturbations. In the other limit, when 
^p- 3> 1, it is the d u 5n s term which dominates, and the Zakharov equation shows density 
rising as ~ t 2 . This is the same scaling as obtained by us using a perturbative approach. The 
density inhomogenities thus created lead to phase mixing and collapse of cavitons. From 
this we infer that a cold 1-D plasma exhibits a 'Langmuir collapse' phenomenon similar to 
what is seen in a warm plasma in 2 or 3 dimensions. The time scale of collapse is of the 
order ~ ( ^4 "^/f^r^ ) 3 plasma periods. Recent experiments on plasma acceleration by laser 
wakefields have shown |J that wave breaking of excited plasma oscillations plays a major 
role in the final acceleration process ; similar physics is likely to be important in the fast 
ignitor concept of laser fusion. We expect the processes discussed in the present paper to 
play some role in such experiments. It should be noted that the simulation results presented 
in this letter are non - relativistic. For many experimental situations the jitter velocity of 
electrons is relativistic and we expect the mass ratio A to be replaced by A e ff ~ m e ///mi 
where m e ff/rrii 3> 1. Under these conditions, the phase mixing effects considered by us 
should become more important. Such investigations are in progress and will be presented 
elsewhere. 
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Figure Captions 

Fig. 1: | 5n d \ n=1 vs. t/T P for A = 1.0 and A = 0.05 
Fig. 2: r mix vs. A and A 
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